Generation of maximum spin entanglement induced by a cavity field 

in quantum-dot systems 
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Equivalent-neighbor interactions of the conduction-band electron spins of quantum dots in the 
model of Imamoglu et al. [Phys. Rev. Lett. 83, 4204 (1999)] are analyzed. An analytical solution 
and its Schmidt decomposition are found and applied to evaluate how much the initially excited dots 
can be entangled with the remaining dots if all of them are initially disentangled. It is demonstrated 
that perfect maximally entangled states (MESs) can only be generated in systems of up to six dots 
with a single dot initially excited. It is also shown that highly entangled states, approximating the 
MESs with good accuracy, can still be generated in systems of odd numbers of dots with almost 
half of them excited. A sudden decrease of entanglement is observed on increasing the total number 
of dots in a system with a fixed number of excitations. 
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I. INTRODUCTION 



Since the seminal papers of Obermayer, Teich, and 
Mahler there has been growing interest in the 
quantum-information properties of quantum dots (QDs) 
in the quest to implement quantum-dot scalable quan- 
tum computers [|[ ^, ^. Those high expectations are 
justified to some extent by recent experimental advances 
in the coherent observation and manipulation of quan- 
tum dots ^ , including spectacular demonstrations of 
the quantum entanglement of excitons in a single dot 
or quantum-dot molecule and observations of Rabi 
oscillations of excitons in single dots Among various 
models of quantum computers based on localized elec- 
tron spins of quantum dots as qubits ||, ^, the scheme 



of Imamoglu ei aZ. |10 is the first, where the interactions 
between the qubits are mediated by a cavity field. This 
approach combines the advantages of long-distance opti- 
cally controlled couplings with long-decoherence times of 
the spin degrees of freedom. Here, we analyze quantum 
entanglement in the Imamoglu et al. model. 

During the last decade, it has been highlighted that 
quantum entanglement, being at the heart of quantum 
mechanics, is also a powerful resource for quantum com- 
munication and quantum-information processing. Quan- 
tum entanglement in interacting systems is a common 
phenomenon. It is obvious that any interacting many- 
body system with defined qubits, if set in a properly 
chosen state, will evolve through states with entangled 
qubits. Surprisingly, quantitative descriptions of the en- 
tanglement dynamics in multiparticle systems are by no 
means satisfactory yet |ll|. Nevertheless, in a special 



case of bipartite entanglement, a number of measures 
have been introduced and studied |l^, . For exam- 
ple, entanglement of a bipartite system in a pure state, 
described by the density matrix pAB = (I^)(^I)aBj can 
be measured by the von Neumann entropy Rq] 

E[pab\ = -Tr{pA log2 pa] = -Tr{/5i3 log2 pb} (1) 

of the reduced density matrix pA — Trslp^s} or, equiv- 
alently, pb — ^t^a{pab}- The entanglement of forma- 
tion of a mixed state of a bipartite system is often mea- 
sured by the so-called concurrence proposed by Hill and 
Wootters[p^. Concurrence has been applied to study en- 
tanglement in various models [ p^ including equivalent- 
neighbor systems [T^. The following two aspects of 
entanglement are especially important: (i) coherent ma- 
nipulation of entanglement and (ii) generation of maxi- 
mum entanglement. The possibility of coherent and se- 
lective control of entanglement in a quantum-dot system 
was analyzed by Imamoglu et al. pO[ |. Here, we would 
like to focus on the latter topic, i.e., the generation of 
the maximally entangled states (MESs) of quantum dots 
in the model of Imamoglu et al. ||l^. MESs are neces- 
sary for the majority of quantum information-processing 
applications. Otherwise, for example, direct application 
of partly entangled states for teleportation will result 
in unfaithful transmission, while superdense coding with 
partly entangled states will cause noise in the resulting 
classical channel. 

The paper is organized as follows. In Sec. II, we de- 
scribe an equivalent-neighbor quantum-dot model and 
give its analytical solution. In Sec. HI, we analyze the 
possibilities of generation of the MESs or their good ap- 
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FIG. 1: Three-level atom in V configuration as a model of a 
semiconductor dot with the conduction-band spin states 
(spin up) of energy sl^'' and |0)n (spin down) of energy si^'' , 
and the effective valence-band state \v)n of energy in the 
nth dot. Key: i^cav, frequency of the common cavity mode; 



lation and creation operators of the cavity mode, re- 
spectively; afn"^ and {ai^^'Y are the corresponding op- 
erators for the laser modes; a^^ is the nth dot operator 

given by af^ = |a;)„„(y|; is the energy of level \x)„ 
{x — 0, l,v); the nth dot levels |0)„ and \v)n are coupled 
by dipole interactions with a strength of g^^ ; analogously, 
g1[^ is the couphng strength between levels |1)„ and |t;)„. 
There is no direct coupling between levels |0)„ and 
in either the same (n = m) or different dots (n ^ m). 
The Hamiltonian (||) simply generalizes, to N dots and 
-I- 1 fields, models of a three-level atom (dot) interact- 
ing with two modes of radiation fields widely discussed in 
the literature (see, e.g., By applying an adiabatic 

elimination method, Imamoglu et al. derived the effec- 
tive interaction Hamiltonian describing the evolution of 
the conduction-band spins of N quantum dots coupled 
by a microcavity field in the form pO| 



H, 



cff 
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„^(A„-A„)t 



frequency of the classical laser field addressed at the nth 



dot; RAt^i^' = Si^^- 
,(1) 



and An = Atj„ — Aa;„ are detunings 



(0) 



cav ; iti-^uJn — L'n 



)*] 



monic mean of AujiP and Auji^^ 



(3) 

in terms of the Pauli spin creation and annihila- 
tion (7^ operators acting on the conduction-band spin 
states of the nth dot. The effective two-dot coupling 
strength between the spins of the nth and mth dots is 
given by Knm{t) = gn(t)gm(t) / where the effective 
single-dot coupling of the nth spin to the cavity field 
{t) I / with Aa;„ being the har- 
For simplicity, the laser 
fields are assumed to be strong and treated classically 
as described by the complex amplitudes El^\t). The 
Hamiltonian \^ was derived by applying adiabatic elim- 
inations of the valence-band states |u)„ and cavity mode 
flcav, which are valid under the assumptions of negligible 
coupling strength, cavity decay rate, and thermal fluctu- 
ations in comparison to ?iA„ and ?iAtj^^^ {x — 0, 1) and 
the energy difference fi^^ — En'^ (see Fig. 1). Moreover, 
the valence-band levels |w)„ were assumed to be far off 
resonance. Although the Hamiltonian (||) describes ap- 
parently direct spin-spin interactions, the real physical 
picture is different: Quantum-dot spins are coupled only 
indirectly via the cavity and laser fields. 

Imamoglu et al. pc| ] applied their model for quan- 
tum computing purposes by implementing the condi- 
tional phase-flip and controUed-NOT (CNOT) operations 
between two arbitrary dots addressed selectively by laser 
fields to satisfy the condition A„ = A^- Here, we are in- 
terested in a realization of an equivalent-neighbor model 
scalable for a large number of dots (even for more than 
100 0). This goal can readily be achieved by assuming 
that all dots are identical and illuminated by a single- 
mode stationary laser field of frequency wi^-* = uj^'"\ 
which implies K,nm{t) = = const. In fact, the condition 
of equivalent-neighbor interactions can also be assured 
for nonidentical dots by adjusting the laser-field frequen- 
interaction Hamiltonian; ficav and a\^^ are the annihi- cies wi^^ to get the same detuning A„ = const, and by 



proximations for different initial conditions of the number 
of excitations and the total number of dots in the system. 



II. QUANTUM-DOT MODEL AND ITS 
SOLUTION 



We will apply the model of Imamoglu et al. |lO[] to de- 
scribe strong equivalent-neighbor couplings of quantum- 
dot spins through a single-mode microcavity field. The 
dots are placed inside a microdisk, put into a microcavity 
tuned to frequency i^cav, and illuminated selectively by 
laser fields of frequencies lo^K Each of N dots with a 
single electron in the conduction band is modeled by a 
three-level atom as shown in Fig. 1. The total Hamil- 
tonian for TV three-level quantum dots interacting with 

-I- 1 quantized fields reads 

H = Hqd + Hp + Hint , (2) 

71 

n 

-Hint = 2^«g„ K + « \ 

71 

n 

where Hqu and Hp arc the free Hamiltonians of the 
quantum dots and the fields, respectively; Hi^t is the 
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FIG. 2: Evolution of the quantum entanglement of E^^{t) 
(solid) and the Schmidt coefficients of P^'^it) (dashed) and 
PI'^ (t) (dot-dashed curves) for systems of A'^ = 2, • ■ ■ , 8 quan- 
tum dots with only one {M — 1) of them initially excited. 
Figure illustrates that the exact maximally entangled states 
can be generated in systems of A'' up to 6 dots only. 



choosing the proper laser intensities \En to obtain the 
effective coupling constants of 5n(i) = const or, equiva- 
lently, Knm(t) = const for every pair of dots. Thus, Eq. 

can be reduced to the effective equivalent-neighbor 
A^-dot Hamiltonian as 



(4) 



where k is the coupling constant. The system described 
by Eq. (^) is sometimes referred to as the spin- 1/2 
van der Waals model |^ , the infinitely coordinated sys- 
tem 1^ , the Lipkin or Lipkin-Meshkov-Glick model , 
or just the equivalent-neighbor model |23| . Let us assume 
that the initial state describing a system of M (M = 
0, ■•■,A^) dots initially excited (i.e., with conduction- 
band spins up) and N — M dots in the ground state 
(conduction-band spins down) is given as 



\m) 



|11_^)a|00_^)b. 

M N-M 



(5) 



Then, we find the solution of the Schrodinger equation 
of motion for the model (Q) in the form 

M' 

m=0 

® {|1)®"|0)^(^^*^^"'}b, (6) 



where M' = min(M, N—M). The states in curly brackets 
{|1)®("-'")|0)®™} denote the sum of all n-dot states with 
{n — m) excitations. For example, {|1)®^|0)®^} stands for 
10011) + 10101) + 10110) + |1001) + |1010) -t- |1100). The 
number of states in the superposition {|l)«'("-™)|o)®™} 
(or equivalently {|1)®™|0)«5("-")}) is equal to the bi- 
nomial coefficient (J^). Thus, for given N and ill, the 
solution (H) contains (^) terms. The energy of the QD 
system described by Eq. (0) is conserved; thus all the 
superposition states in Eq. ^) have the same number M 
of excitations. We find the time-dependent superposition 
coefficients in Eq. (^) as 

M' 



Cl"{t) - exp{zKA^ + l-n) 

- M{N - M)]nt) (7) 



n=0 



in terms of 



uNM 



fc=0 



N -2k 
M -k 



N + l-2k\ J N -2k 
n — k I \n — k — 1 



(8) 



where are binomial coefficients. Our solution can 
be represented in a biorthogonal form via the Schmidt 
decomposition 



M' 



(9) 



m=0 



where |(/)Tn(i))A and \ipm{t))B are the orthonormal basis 
states of subsystems A and B, respectively. We find that 
the real and positive Schmidt coefficients can be related 
to the squared module of superposition coefficients as 
follows: 



NM 



it) 



M\ fN - M 



c^'m^ (10) 



while the phases of C^^^^{t) are absorbed into the def- 
inition of the basis states \4'm{t)) A and \^m{i))B- The 
Schmidt coefficients are normalized to unity. The evo- 
lutions of all P^^^ for systems with single and two ex- 
citations are given in Figs. 2 and 3, respectively. We 
observe that the evolution of Schmidt coefficients is pe- 
riodic with the period of kT = 2tt/N for systems with a 
single {M — 1 or, equivalently, M = N — 1) excitation 
(Fig. 2), and tt periodical {2tt periodical) for systems of 
even (odd) numbers of dots with higher numbers of exci- 
tations (see Fig. 3). For brevity, only half of the period 
is depicted in the right-hand panels of Fig. 3. 



III. ENTANGLEMENT IN QUANTUM-DOT 
SYSTEMS 



We address the following questions: How much can 
the initially excited dots (say, subsystem A) be entangled 



3 





1.6 


Ie1.2 






> 


0.8 




0.4 


w 






0.0, 




1.6 


Is 1.2 


Ph 




> 


0.8 


S 
z 


0.4 


U 






0.0, 




1.6 


Is 1.2 


Ph 






0.8 


NM 


0.4 








0.0, 



N=4 M=2 



N=5 M=2 



N= 10 




N=20 















'■\ 1 / 












■■■■ \ ■*■ 











n 

Kt 



3n/2 



FIG. 3: Evolution of the entanglement of S^^(t) (solid) 
and all Schmidt coefficients Pa^{t) (dashed), Px{t) (dot- 
dashed), and P2^^(t) (dotted curves), in systems with two 
(M = 2) dots initially excited. 
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FIG. 4: Maximum entanglement E„ 



maxt£;'^"(i) (solid 



bars) , measured in ebits, as a function of the excitation num- 
ber M generated in systems of = 10, 20, 30 and 31 dots. 
The empty staircase corresponds to entanglement of i?MES 
for the MESs. The figure illustrates that the highest entan- 
glement, closest to S^gg, can be generated in systems with 
M = [A^/2] excitations. On decreasing M or (iV - M), the 
entanglement decreases. The discrepancy between E^^^ and 
^MES becomes more pronounced with increasing A'^ especially 
for < M < [N/2]. 



with the remaining dots (subsystem B) in the equivalent- 
neighbor system of initially all disentangled dots if the 
evolution is governed by Hamiltonian (^)? And whether 
the maximally entangled states can be generated exactly 
or, at least, approximately in systems of an arbitrary 
number N of dots while M of them are excited. 

With the help of an explicit form of the Schmidt de- 
composition, it is convenient to calculate the entangle- 
ment (|l|) via the Shannon entropy 

E^''^{t) ^ Emt)){m\] 

M' 

= -E^™*'wi°g2nrw (11) 

m=0 



of the Schmidt coefficients given for our system by (10). 
By applying Eq. (jl^), we can determine the maxi- 
mum entanglement given by E^^(t) = maxj i?^*^(i), 
which can periodically be generated during the evolu- 
tion of TV-dot system with M excitations. The coef- 
ficients (H^), as well as (0), possess the symmetry of 
P^*^(i) — P,^'^~*^(t), which impHes equal evolutions 
of entanglement 



{12) 



in the N-dot systems with M and N — M excitations. 
Figure 4 shows this symmetry in a special case for maxi- 
mum entanglement of max* E'^^''{t) = maxt E^'^-^{t). 

To solve the second problem proposed at the begin- 
ning of this section, we have to determine the quantum 



correlations of the maximally entangled state of two sub- 
systems having d equally weighted terms in its Schmidt 
decomposition. According to the theorem of Bennett et 
al. [Q, the MES has log2 d ebits of entanglement, where 
d is the Hilbert space dimension of the smaller subsys- 
tem. Thus, in our case, the MES of the subsystem A 
consisting of M dots and the subsystem B oi N — M 
dots has 



<ES = log2[niin(M,iV-M)-M] 



(13) 



ebits of entanglement. In particular, the MES in the TV- 
dot system with a single initial excitation has only 1 ebit 
independent of N . The empty staircase in Fig. 4 and 
solid lines in Fig. 5 correspond to i^MES- ^'-^ show a 
deviation of a given state from the MES, it is convenient 
to use the relative (or scaled) entanglement defined to be 



NM 



E, 



NM 



T^NM 
^MES 



max 

t E, 



E^^'\t) 



NM 
MES 



In the simplest nontrivial case, for M 
coefficients reduce to 



(14) 

1, the Schmidt 



i(t) =4 



7V2 



N 

-Kt 



(15) 



and P^^{t) 



1 — Pi (t), which enable a direct calcu- 
lation of the entanglement E^^{t) with the help of Eq. 
(|rT|). The evolutions of entanglement and the Schmidt 
coefficients of Pj^^{t) for m=0,l, are depicted in Fig. 2. 
The quantum-dot systems evolve into the MESs at evo- 
lution times that are the roots of the equation 
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FIG. 5: Maximum entanglement -Eniax as a function of the 
total number A'^ of dots generated in systems with M=l,2,3, 
and excitations. The solid lines and empty staircase 

correspond to -Emes- On the scale of the figure, an apparent 
plateau occurs for smaller than some critical value Nm . For 
higher than Nm and fixed M, a monotonic decrease of the 
maximum entanglement is clearly visible. One concludes that 
arbitrary high entanglement can be achieved by increasing A'' 
and keeping half M = [A'^/2] of the system excited. 
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FIG. 6: The same as in Fig. 5 but for the relative maximum 
entanglement e^f^ = ^-majf /^'mes • The figure shows that 
the apparent plateau for finite M actually occurs for M = 1 
only. The first and second highest maxima of entanglement 
correspond to A'' equal to 2M + 1 and 2AI + 5, respectively. 



= £:^A'^(i) = 2K.^^^sm{NKt) 



X log 



2 fN 

— 7 CSC Kt 



- 1 



Thus, we get 



^i'^larccscf ^v/2(iV^l) 



(.16) 



(17) 



and Kt" — tt/N. We find that the maximum entangle- 
ment, equal to E^^{t') — 1 ebit, can be achieved at evo- 
lution times t' for < 6 only. For A^ > 6, a real solution 
for t' does not exist. Another explanation of this result, 
as illustrated in Fig. 2, can be given as follows: The 



maximum entanglement corresponds to the Schmidt co- 
efficients mutually equal or, in general, the least different. 
But the MES corresponds solely to the former case. As 
seen in Fig. 2, the condition P(f ^(^') = Pf^it') is strictly 
satisfied for A^ < 6. The entanglement for A^ > 6 reaches 
its maximum at evolution times t" . This maximum value 
is given by 



-2(Af-l)log2[4(A^-l)]}, 



(18) 



which is less than unity and monotonically decreases with 
increasing A^ as clearly illustrated in Figs. 5 and 6 for 
M=l. Thus, the perfect MESs cannot be generated in 
systems of A^ > 6 dots. Nevertheless, a good approxima- 
tion of the MESs can also be obtained for N = 7. On 
the scale of Fig. 2, maxt E'^ (t) = E'^'^{t:/7) ^ 0.9997 
is close to unity since P^^{Tr/7) and P-l^^ {Tr/7) are al- 
most the same. It is worth noting that a critical value of 
A^ = 6 was also found, although in the different context 
of the pairwise entanglement measured by the concur- 
rence 14 , for an equivalent-neighbor model of entangled 
webs in Ref. [ p6| . In comparison, a critical value of 
A^ = 6 for the concurrence in the equivalent-neighbor 
isotropic or anisotropic Heisenberg models was not ob- 
served (see, e.g., [|l^). Similarly, generation of the MESs 
in an equivalent-neighbor quantum-dot model of Reina 
et al. was discussed only in two special cases of the Bell 
{N=2) and Greenberger-Horne-Zeifinger (GHZ) {N=3) 
entangled states [|l^. Thus, no critical behavior of en- 
tanglement as a function of A'^ was reported there. 

The case for M = 1 is the only one where the general 
formula ( [To| ) for the Schmidt coefficients simplifies to a 
compact form for arbitrary evolution times. Thus, for 
clarity, we present mainly numerical results for M > 2. 
For example. Fig. 3 illustrates that the exact MESs can- 
not be generated in systems with M — 2 excitations at 
any evolution time. This conclusion can be drawn from 
the observation that P,^^(t) for to=0,1,2 do not cross si- 
multaneously at any times in the period. Nevertheless, 
the MESs can be approximated with good precision. The 
highest possible entanglement, corresponding to the least 
mutually different Pj^'^, is observed for A^ = 5 and 9, 
where the relative entanglement deviates from unity at 
the order of 10~^ and 10"'*, respectively (see Fig. 6 for 
M=2). The states generated in A^-dot systems with three 
excitations can be entangled up to e^^^ = 0.9996 (first) 
and ej^^x = 0.9990 (second maximum) for the relative 
entanglement (see Fig. 6 for M=3). It is interesting to 
compare the relative entanglement of e^^^, depicted in 
Fig. 6, with the "absolute" entanglement of E^^^ pre- 
sented in Fig. 5. By analyzing the numerical data given, 
in part, in Fig. 6, we find the following rule: The maxi- 
mally or almost maximally entangled states can be gener- 
ated in systems of A^ = 2M + 1 dots with M excitations. 
Slightly worse entanglement can be achieved in systems 
of A^ = 2AI + 5 dots with Af excitations. Thus, systems 
composed of odd rather than even numbers of dots enable 
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FIG. 7: The inverse of the maximum entanglement, {E^ 
(dots) measured in ebits"^, and its approximation (solid lines) 
as a function oi N > 2M + 5 generated in systems with M — 
1, 2, 3 excitations. 



generation of the entangled states better approximating 
the MESs for M > 1. This is clearly illustrated in Fig. 6 
for M = [N/2], i.e., the integer part of N/2. We observe 
that the system of odd and large numbers {N > 2M + 5 
for M > 1) of dots is the most entangled at evolution 
times Kt = {1 + 2k)TT for fc = 0, 1, ... (see, e.g.. Fig. 3 for 
iV=ll). In this special case, the Schmidt coefficients can 
be written compactly via 



C, 



NM 



= 2"m!(7V-2Af) 



(A^-2to-2)! 

m. 



(19) 



For nt — kn and even N, in contrast to odd N, the 
entanglement vanishes. The maximum entanglement of 
E^^^ for TV > Nm = 2M + 5- Sim can be well fitted by 
the inverse of linear functions as shown in Fig. 7. 



IV. CONCLUSION 

We studied the evolution of the conduction-band spins 
of quantum dots in the model of Imamoglu et al. mn . We 



found the analytical solution and its Schmidt decomposi- 
tion for the equivalent-neighbor model and applied them 
in our study of bipartite entanglement in quantum-dot 
systems with arbitrary numbers of dots and their exci- 
tations. We have raised and solved the problem to what 
extent the initially excited dots can be entangled with 
the remaining dots if all of them are initially disentangled 
in the equivalent-neighbor energy-conserving model. We 
have shown that the perfect maximally entangled states 
can only be generated in systems of N = 2, • ■ • , 6 dots 
with a single dot initially excited. Nevertheless, highly 
entangled states, being excellent approximations of the 
MESs, can periodically be generated in systems of odd 
numbers N of dots with the number M of excitations 
equal to M = {N — l)/2 (leading to the best approxima- 
tion) and M = (A^ — 5)/2 (giving a slightly worse approx- 
imation). If we increase N beyond Nm = 2M -f 5 — Sim, 
the entanglement decreases monotonically as described 
by the inverse of linear functions. 



ACKNOWLEDGMENTS 



We thank J. Bajer, T. Cheon, T. Kobayashi, H. Mat- 
sueda, and I. Tsutsui for their stimulating discussions. 
Y.L. acknowledges support from the Japan Society for 
the Promotion of Science (JSPS). This work was sup- 
ported by a Grant-in-Aid for Scientific Research (B) 
(Grant No. 12440111) and a Grant-in-Aid for Encour- 
agement of Young Scientists (Grant No. 12740243) by 
the Japan Society for the Promotion of Science. 



[1] K. Obermayer, W. G. Teich, and G. Mahler, Phys. Rev. 
B 37, 8096 (1988); 37, 8111 (1988). 

[2] A. Barenco et al, Phys. Rev. Lett. 74, 4083 (1995); 
H. Matsueda and S. Takeno, lEICE Trans. Fundamen- 
tals E79-A, 1707 (1996); H. Matsueda, Superlattices Mi- 
crostr. 24, 423 (1998); G. D. Sanders et al., Phys. Rev. 
A 60, 4146 (1999); E. Biolatti et al, Phys. Rev. Lett. 85, 
5647 (2000). 

[3] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998); G. Burkard et al, Phys. Rev. B 59, 2070 (1999). 

[4] S. N. Molotkov, JETP Lett. 64, 237 (1996); S. Bandy- 
opadhyay, Phys. Rev. B 61, 13 813 (2000); T. Ohshima, 
Phys. Rev. A 62, 062316 (2000); J. Levy, ibid. 64, 052306 
(2001). 

[5] N. H. Bonadeo et al, Science 282, 1473 (1998). 



[6] T. H. Oosterkamp et al.. Nature (London) 395, 873 
(1998). 

[7] G. Chen et al. Science 289, 1966 (2000). 
[8] M. Bayer et al.. Science 291, 451 (2001). 
[9] T. H. Stievater et al, Phys. Rev. Lett. 87, 133603 (2001); 
H. Kamada et al., ibid. 87, 246401 (2001). 
[10] A. Imamoglu et al, Phys. Rev. Lett. 83, 4204 (1999); 

Fortschr. Phys. 48, 987 (2000). 
[11] M. Murao et al, Phys. Rev. A 57, R4075 (1998); A. 
Thapliyal, ibid. 59, 3336 (1999); C. H. Bennett et al, 
ibid. 63, 012307 (2000); W. Diir et al, Phys. Rev. Lett. 
83, 3562 (1999); N. Linden et al, ibid. 83, 243 (1999); 
H. A. Carteret et al. Found. Phys. 29, 527 (1999). 
[12] C. H. Bennett et al, Phys. Rev. A 53, 2046 (1996). 
[13] V. Vedral et al, Phys. Rev. Lett. 78, 2275 (1997); V. 



6 



Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 (1998). 
[14] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 

(1997); W. K. Wootters, ibi d 80. 2245 (1998). 
[15] W. K. Wootters, e-print |quant-ph/00011l4 K. M. 

O'Connor and W. K. Wootters, Phys. Rov. A 63, 052302 

(2001); D. Gunlycke et al, ibid. 64, 042302 (2001); M. 

C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 87, 

017901 (2001). 

[16] M. Koashi, V. Buzek, and N. Imoto, Phys. Rev. A 62, 
050302 (2000). 

[17] X. Wang and K. M0lmer, Eur. Phys. J. D 18, 385 (2002). 
[18] J. H. Reina et al, Phys. Rev. A 62, 012305 (2000). 



[19] 



[20] 

[21] 
[22] 

[23] 



C. C. Gerry and J. H. Eberly, Phys. Rev. A 42, 6805 
(1990); M. Alexanian and M. Bose, ibid. 52, 2218 (1995); 
Y. Wu, iUd. 54, 1586 (1996). 

R. Dekeyser and M. H. Lee, Phys. Rev. B 19, 265 (1979); 

43, 8123 (1991); 43, 8131 (1991). 

R. Botet et al, Phys. Rev. Lett. 49, 478 (1982). 

H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 

62, 188 (1965). 

J.-M. Liu and C. Miiller, Phys. Rev. A 42, 5854 (1990); 
Phys. Rev. B 44, 12 020 (1991). 



published in the Physical Review A 65 (2002) 062321 
and selected to 

Virtual J. Nanoscale Sci. Tech. ( http://www.vjnano.org/nano/ ) 6 (2002) Issue 1; 
Virtual J. Quantum Information ( ht tp : // www . v j quantuminf o . or^ ) 2 (2002) Issue 7. 



7 



